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Exciton dispersion as a function of center-of-mass momentum, Q, is essential to the understanding 
of exciton dynamics. We use the ab-initio GW-Bethe Salpeter equation method to calculate the 
dispersion of excitons in monolayer M 0 S 2 and find a nonanalytic light-like dispersion. This behavior 
arises from an unusual \Q\ term in both the intra- and inter-valley exchange of the electron-hole 
interaction, which concurrently gives rise to a valley quantum phase of winding number two. A 
simple effective Hamiltonian to order with analytic solutions is derived to describe quantitatively 
these behaviors. 


PACS numbers: 73.22.-f, 71.35.-y, 78.67.-n 

An exciton is a neutral excitation of an interacting 
electron system, consisting of a bound excited electron 
and hole pair with energy depending on its center-of- 
mass momentum Q. In crystals, the exciton dispersion 
relation (energy vs. Q) forms an exciton band struc¬ 
ture. In optical absorption, the momentum of a photon 
is converted to that of an exciton, resulting in a small Q 
exciton for visible light. Recent model Hamiltonian cal¬ 
culations predict that the dispersion of the lowest energy 
optically active exciton bands in monolayer M 0 S 2 form 
a Dirac cone (as in graphene) due to the exchange inter¬ 
action coupling electron-hole states in different valleys 
(the intervalley exchange) [1]. Another study based on a 
tight-binding formulation of the BSE approach suggested 
also linear dispersion for the lowest energy excitons near 
Q = 0, although the dispersion relation (on a 45 x 45 
finite-Q grid) was not resolved at the lengthscale of the 
momentum of light [2]. The effects of the interplay of 
inter- and intra-valley exchange and local fields are not 
explicitly investigated in these model calculations, which 
also miss terms in both Q and orders. 

In this work, we calculate the exciton dispersion of 
M 0 S 2 , a prototypical transition metal dichalcogenide 
(TMD), from first principles, using the ab-initio GW- 
Bethe Salpeter equation (GW-BSE) method [3l[4], and 
find a highly unusual low-energy dispersion consisting of 
a nonanalytic v-shaped upper band that is degenerate 
with a parabolic lower band at Q = 0, consistent quali¬ 
tatively with previous tight-binding results [2] . We show 
that the physical origin of this highly nonanalytic be¬ 
havior (the |(5|-dependence) comes from |(5|-dependent 
terms in both the intravalley and intervalley exchange 
interaction, which arise from the unique electronic struc¬ 
ture and the quasi-2D nature of the Coulomb interaction 
in atomically thin TMDs. Local-field effects introduce 
additional interaction terms and are responsible for the 
splitting of optically bright and dark excitons. Moreover, 
the theory gives a valley quantum phase of winding num¬ 
ber two (or chirality two), which we show should manifest 
in optical experiments as a phase difference between the 


longitudinal and transverse response. A similar winding 
number is found in Ref. [1] , although in this model Hamil¬ 
tonian study, intravalley exchange is neglected resulting 
in a Dirac cone dispersion. 

Transition metal dichalcogenides are layered, weakly- 
coupled materials. In the 2H monolayer form, TMDs 
are direct band gap semiconductors with the gap at 
the K and points of the 2D hexagonal Brillouin 
zone Eli- Due to a lack of inversion symmetry and pres¬ 
ence of three-fold rotational symmetry, a valley selective 
response emerges: the K and valley couples respec¬ 
tively only to left or right circularly polarized light IZHIQ]. 
Moreover, as a result of confinement and reduced screen¬ 
ing, excitons have very large binding energies mHni. 

The quasiparticle (QP) bandstructure of monolayer 
M 0 S 2 calculated within the GW approximation O [15] 
(see Supplemental Materials) is shown in Fig. The 
direct gap of 2.67 eV is at the K and points. Spin- 
orbit coupling splits the valence band edge by 147 meV 
and the conduction band edge by 3 meV. Over an ex¬ 
tended range in the K and valleys, the spin of the 
electron along the direction perpendicular to the layer, is 
a good quantum number [8|. Thus, the concept of singlet 
and triplet exciton states is no longer well-defined. In¬ 
stead, the electron-hole Hamiltonian or BSE matrix can 
be decoupled into transitions between bands of like-spin 
(Fig. [^b-d)) that are optically allowed and transitions 
between bands of unlike spin that are optically forbid¬ 
den. We calculate the exciton dispersion of both the like- 
spin transition states and the unlike-spin transition states 
within the BSE formalism following Ref. [4]. The result¬ 
ing exciton dispersion near Q=0 and Q=K is shown in 
Fig-i Fig. only presents exciton states involving pre¬ 
dominantly transitions from the topmost valence band 
(the V f K and v K' bands in Fig. [^, which is the 
so-called A series in the literature. Similar results are 
obtained for the spin-orbit split B series (see Supplemen¬ 
tal Materials). 

In Fig. near Q = 0, the lowest energy exciton com¬ 
plex shown is the ”ls”-like states of the A series excitons. 
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FIG. 1. (Color online) (a) QP bandstructure of monolayer M 0 S 2 . Lowest energy transitions corresponding to electron-hole 
excitations with momentum transfer of Q = 0 and Q = X are shown with the labeled arrows, (b-c) Schematics of interactions 
between two electron-hole pairs with momentum near Q = 0, corresponding to BSE matrix elements {vckCl\H^^^\v'c''k.'Cl) 
for: (b) two like-spin transitions within one valley and (c) two like-spin transitions from different valleys, (d) Schematic of two 
electron-hole pairs with momentum near Q = iF giving rise to BSE matrix elements for like-spin transitions. 


which has a binding energy of 0.63 eV for the like-spin 
transition states and 0.65 eV for the unlike-spin transi¬ 
tion states. At Q = 0, the two-fold degeneracy of the 
exciton states due to K and iF’ valley degeneracy is pro¬ 
tected by the time-reversal symmetry. Away from Q = 0, 
for optically active like-spin transition states, the degen¬ 
erate bands split as |(5| increases. Consistent with previ¬ 
ous model calculations [U [2j |27l [28], this splitting of the 
like-spin transition states is due to intervalley exchange. 
The exciton bands for unlike-spin states, which have zero 
exchange interaction, remain doubly degenerate (the blue 
lines in Fig.[^. However, contrary to previous model pre¬ 
dictions [ 1 ] , our ab-initio results do not find a Dirac cone 
near Q = 0 in the dispersion of the like-spin transition 
”ls” bands. In fact, we find a ”v-shaped” nonanalytic up¬ 
per band and a parabolic lower band, with both bands 
increasing monotonically with \Q\ down to the smallest 
sampled \Q\ of 4 x 10“^which corresponds to 0.3% 
of the distance from E to iF in the Brillouin zone. 


We also calculate the dispersion of the intervalley ex- 
citons (Fig. I^b)) —i.e., an exciton with the electron and 
hole in different valleys (with Q = K )—and find that 
the unlike-spin transitions state is slightly lower in en¬ 
ergy than the similar exciton at Q = 0 as a consequence 
of the spin-orbit splitting of the conduction bands and the 
screened Coulomb interaction (see Supplemental Materi¬ 
als). 


We now provide the physical origin of the nonana¬ 
lytic behavior of the dispersion seen in Fig.j^a) (closeup 
in Fig. 1^). The exciton dispersion is obtained from 
first principles by solving the Bethe-Salpeter equation 
(BSE) HES] for electron-hole pair excitations with fi- 
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EIG. 2. (Golor online) Exciton dispersion of monolayer M 0 S 2 
near (a) Q=0 and (b) Q=K along the K — T — K' direction. 
Red (blue) lines indicate states arising from like (unlike)-spin 
transitions. The label A refers to states involving transitions 
from the highest valence band at K/K’. ”B”-states involving 
transitions from the second highest valence band are similar 
and shown in the Supplemental Materials. 


nite center-of-mass momentum Q 

(-£^ck+Q ~ -^vk)A^ckQ + 

v'c'k' 

Here, S indexes the exciton states; is the ampli¬ 

tude of the free electron-hole pair consisting of an elec¬ 
tron in |ck -f Q) and one missing from |'ck); Dq is the 
exciton excitation energy; F^ck+Q and Ey\^ are the quasi¬ 
particle energies, and is the electron-hole interac¬ 
tion kernel. The kernel consists of a direct term and an 
exchange term[4]: {vckQ\K^^\v'EliQ) = {vckQ\K^ + 
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K^\v 'The exchange term is 
(i;ckQ|i^"^|i;'c'k'Q) = 

^Mc„(k, Q, G)^;(Q + G)M,%,(k', Q, G), (2) 

G 

where G are reciprocal lattice vectors; v is the 
bare Coulomb interaction, and M is defined as 
M^n/(k, Q, G) = (nk + Q|e*^^+^^'^|n'k). The exchange 
term is only non-zero for like-spin transitions (i.e. exci- 
tons with total spin along z equal to zero). The direct 
term is 

{vckQ\K'^\v'c'k'Q) = 

- M*,,(k + Q,q,G)W^GG'(q)M„„.(k,q,G), (3) 

GG' 


where q = k — k^ and W is the screened Coulomb inter¬ 
action. 

The solutions of the BSE matrix are exciton states 
I^q) = S^;ck^?ckQl'^<^kQ)- subspace spanned 

by the two lowest energy states, at exactly Q = 0, the 
exchange matrix elements are also diagonal. At Q = 0, 
the exchange term is the diagonal constant 

C ^ (Soms'o) = 5ss'27re^ Y1 E 

Gt^O vckv'c'k' . . 

^vckO^v'c'k'O 1^1 ’ 


where for small Q’s in 2D, v{Q + G) = ~ IW' 

This constant term is 20 meV in our ab-initio calcula¬ 
tion. It is purely a local-field effect and is responsible 
for the splitting between the like-spin(bright) and unlike- 
spin(dark) states at Q = 0. 

We now derive an effective Hamiltonian to 

describe the main physics for the Is complex in a basis of 
”excitonic” functions from the individual valleys in the 
tight-binding limit, given in Eq. below. At Q = 0, 
both the like and unlike-spin excitonic levels are doubly 
degenerate, with the amplitude of one state {\Sq)) 

confined to the K valley and that of the other {\Sq )) 
confined to the valley. Near Q = 0, it is sufficient to 
use the following basis functions (which are of the Bloch 
form of excitons from a specific valley in the tight-binding 
limit) to expand the true exciton state |5'q): 










where R = , {o)q and (i)q are pseudospinors de- 

noting respectively basis functions on K and K' valleys, 
and = E,,cfc ^Ikol^ckQ). 

In this basis, is a 2 x 2 matrix, and the in¬ 

travalley exchange term (matrix element between basis 
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FIG. 3. (Color online) (a) Close-up of dispersion of Ais near 
Q=0. Ab-initio values are stars. Fit to effective Hamilto¬ 
nian (Eq. pr] ) are solid lines. Red(blue) lines indicate states 
arising from like (unlike)-spin transitions, (b) Valley pseu¬ 
dospin texture of the (upper) nonanalytic like-spin transition 
band around Q=0 for states of fixed energy in Q-space. (c) 
Valley pseudospin texture of the (lower) parabolic like-spin 
transition band, (d) Optical absorbance of linearly polarized 
light at fixed incidence as the polarization vector e is rotated 
over 360°. The angle of the polarization vector, 0', is defined 
with respect to the vector formed by the intersection of the 
polarization plane (blue) and the xy-plane. Red (black) indi¬ 
cates the absorbance of states arising from the lower (upper) 
like-spin band, (e) Energy difference between the upper and 
lower like-spin bands that is probed as 0, the angle between 
wavevector of light and the z-axis, is changed, for light of 
huj ~ 2 eV. The inset shows how 0, 0, and 0' are defined. 


functions in the same valley) is 


{S^\K^K) 


E ^vcko^v^c'k'o (vckQ Ilu'c'k'Q). 

vckv'c'k' 

(6) 


Using a Q - p expansion of the QP states in the M matrix 
elements (Eq. to second order in Q • p, the intravalley 
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exchange term (Eq. becomes 
{S^\K-\S^) = {S^'\K-\S^')= Y1 ^Iko^^vk'o 

vckv'c'k' 

{Q • a((u„k|p|wck))*Q • a'(w«'k' |p|Mc'k')t’(Q) 

( 7 ) 

Here, we have separated contributions from the G = 0 
and G 7 ^ 0 Fourier components (Higher order terms are 
described explicitly in the Supplemental Materials). We 
note that a, a', and c are factors that depend on the QP 
states but are independent of Q [30] . In Eq. the G = 0 
term, has a nonanalytic dependence on |(5| 

IQ-E“^«cko(w«k|p|wck)|^f(Q) oc^lQI- (g) 

vck 

Here, H is a proportionality constant, and we have made 


use of the fact that for small Q in 2D, '^(Q) 


27re^ 

IQI ’ 


and the fact that (0|p|5'^) oc x ^ iy due to the Cs 
symmetry [7], which eliminates the dependence on the ori¬ 
entation of Q. Up to order the intravalley exchange 
(Eq.§ has the form 


{S^\K^\S^) = {S^'\K-\S^') = C + A\Q\+(3Q^. (9) 


C is the splitting of bright and dark states given in Eq. 
and is a proportionality constant that is a real number, 
which arises from the local fields. Likewise, the interval¬ 
ley exchange may be shown to be 


{S^\K^\S^') = {{S^'\K^\S^)r = A\Q\e-^^o + I3'Q\ 

( 10 ) 

where 0 is the angle of Q defined with respect to the 
x-axis and is a complex number. 

After performing a similar analysis of the direct term, 
we find an effective Hamiltonian 


H 


BSE 


h 

(Q) = Uoff + A[1 + cos{20)(Jx + sm{20)(jy]\Q\ + [( 2 ^ + + /3)1 + P\cfx + 


( 11 ) 


where Qq is the excitation energy of the exciton with 
Q = 0; M = rrie ^ rrih is the QP band mass of the free 
electron-hole pair at LC or LC'; a is a constant from the 
order of contribution from the direct term; /3 and 
are constants from the order of contribution from the 
intravalley and intervalley exchange respectively; (Tx and 
Gy are the Pauli matrices (see Supplemental Materials for 
detailed derivation). An effective Hamiltonian may be 
obtained for the other exciton complexes through similar 
analyses. 


We fit Eqs. 12 and 13 to our ab-initio calculation 


Diagonalizing the effective Hamiltonian (Eq. 11) gives 
us two solutions: one with a parabolic dispersion 


U_(Q) = Oo 


+ (^ 2 M + ^ ^ ( 1 ^) 


and the other with nonanalytic dispersion 

^+(Q) = Uo + 2A\Q\ + + a + ^ + |/3'|^ (13) 

Hence, it is the combination of intervalley and intravalley 
exchange that results in a lower band with a parabolic 
dispersion (Eqjl^ and an upper band with a nonanalytic 
dependence on jQI (Eq.[^, as seen in our ab-initio calcu¬ 
lation (Fig. [^). We emphasize that this |(3|-dependence 
at small Q is a consequence of the 2D Coulomb inter¬ 
action and the opposite angular polarization of the elec¬ 
tronic states in the K and K' valleys, which removes the 
dependence on the direction of Q. 


(Fig. I^a)) to obtain the values of the proportionality 
constants. A, a, (3, and /3'. We find A = 0.9 eV-A. 
The dominant term comes from the QP effective mass 
. M for monolayer M 0 S 2 is roughly l.lmo- From 
the dispersion of the unlike-spin states, we find that 
a = —0.9 eV-A^, and the effective mass of the unlike- 
spin transition(dark). Is exciton is roughly M* = 1.5mo. 
The effective mass of the parabolically-dispersing like- 
spin(bright). Is exciton is M* = 1.4mo, indicating that 
(3 — \P'\, the difference between quadratic terms in the 
intervalley and intravalley exchange, is small, only about 
0.2 eV-A^, while the magnitude of /3 is about 4 eV-A^. 
Our results demonstrate an enhancement of about 30% 
in the center-of-mass effective mass M* of the parabolic 
band excitons due to electron-hole interaction m- 
Since the difference between the magnitude of the 
quadratic terms in the intervalley and intravalley ex¬ 
change, P — \P'\, is small, we may to good approxima¬ 
tion neglect this difference (i.e. we take P' = 

See Supplemental Materials.). Then, the solution for the 
parabolic like-spin transition exciton band simplifies to 




-ie\QK 




(14) 


and the solution for the nonanalytic like-spin transition 
exciton band simplifies to 






Q/ 


\S^ ))• 


(15) 








5 


For both states, the valley pseudospin winding number 
around Q=0 is 2. The pseudospin texture of these states 
in a circle in Q-space is shown in Fig.[^b-c). To study the 
optical response, we project the momentum and polar¬ 
ization of incident linearly-polarized light into the layer 
2D plane. The momentum Q of the excited exciton is 
equivalent to the in-plane component of the wavevector 
of the light. The in-plane projection of the electric field 
polarization vector can be decomposed into components 
perpendicular (transverse) to and parallel (longitudinal) 
to the momentum transfer Q. We define the angle of the 
polarization vector with respect to the transverse pro¬ 
jection as 0'. The absorbance of the upper and lower 
bands with respect to 0' has a phase difference of 90°, 
and the intensity of each band peaks twice as the po¬ 
larization angle is rotated over 360° (Fig. and Sup¬ 
plementary Materials). The different pseudospin texture 
results in the upper band coupling to the longitudinal 
component of the in-plane projection of the electric field 
and the lower band coupling to the transverse component 
of the in-plane electric field (see Supplemental Materi¬ 
als). Since the transverse projection is always in the 2D 
plane, only the intensity of optical absorbance from the 
upper band changes with the angle of the incident light 
with respect to the z-axis (0). Intensity is maximum at 
normal incidence (0 = 0) and minimum at grazing in¬ 
cidence (0 ~ 90°). The energy difference between the 
upper and lower band for different Q on the dispersion 
curve can also be probed by changing 0. For photons 
with huj ^ 2 eV, the largest energy difference is about 
1.5 meV when the wavevector of the light is nearly par¬ 
allel to the plane. 

In summary, we have computed the exciton dispersion 
of M 0 S 2 from first principles and find an unusual disper¬ 
sion with a parabolic lower band and a v-shaped upper 
band for the lowest energy like-spin transition exciton 
complex near Q = 0. This dispersion is due to the in¬ 
terplay of the intervalley and intravalley exchange, both 
of which have a |(5|-dependent behavior. We have de¬ 
rived a simple effective Hamiltonian and analytic solu¬ 
tions describing this physics and predict that the split¬ 
ting of the exciton bands can be measured with a lin¬ 
early polarized optical beam tilted away from normal in¬ 
cidence. We expect any 2D semiconductor with excitons 
with amplitude concentrated in a small portion of the 
Brillouin zone to exhibit similar nonanalytic exciton dis¬ 
persion near Q = 0. We also show that interaction effects 
increase the exciton mass M* by > 30%. First-principles 
results for the intervalley excitons with Q = if and other 
finite Q excitons are also obtained. 
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